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Measurements of the linear power spectrum of galaxies have placed tight constraints on neutrino 
masses. We extend the framework of the halo model of cosmological nonlinear matter clustering 
to include the effect of massive neutrino infall into cold dark matter (CDM) halos. The magnitude 
of the effect of neutrino clustering for three degenerate mass neutrinos with m^^ = 0.9 eV is of 
order within the potential sensitivity of upcoming weak lensing surveys. In order to use these 

measurements to further constrain — or eventually detect — neutrino masses, accurate theoretical pre- 
dictions of the nonlinear power spectrum in the presence of massive neutrinos will be needed, likely 
only possible through high-resolution multiple particle (neutrino, CDM and baryon) simulations. 

PACS numbers: 98.80.-k,98.65.-r,14.60.Pq 



I. INTRODUCTION 



The detection of neutrino flavor oscillation in conjunc- 
tion with cosmological arguments has highly constrained 
neutrino mass eigenvalues. Solar neutrinos Q, Ql and 
atmospheric neutrinos 0, Q oscillate from one flavor to 
another. The KamLAND reactor neutrino detector has 
found evidence for neutrino oscillations consistent with 
the inferred solar neutrino oscillation parameters Q . The 
K2K long-baseline experiment has found evidence for 
neutrino oscillations consistent with the atmospheric re- 
sults While flavor oscillation experiments constrain 
only the neutrino mass differences, cosmological argu- 
ments have the advantage of constraining the total mass. 
The present cosmological upper limits are competitive 
with terrestrial experiments and expected to improve 
substantially with time. 

Massive neutrinos influence the large scale structures 
of the universe in a well-deflned way (3] because they do 
not cluster, thereby reducing the amount of matter that 
can accrete into potential wells. The galaxypower spec- 
trum has been measured on large scales 0, leading to 
upper limits on the sum of neutrino mass ranging from 
0.7 to 1.8 eV [13, HI d El (depending on assumptions 
and data sets). Combining estimates of the linear mat- 
ter power spectrum from the Lyman-a forest in the Sloan 
Digital Sky Survey (SDSS) with estimates of the bias of 
galaxies in the SDSS with galaxy-galaxy lensing, a tight 
limit of 0.42 eV is inferred for the 95% C.L. limit on the 
sum of three degnerate mass neutrinos |14| . The approx- 
imation of a linear spectrum that is valid for such large 
scale measurements is no longer valid on smaller scales 



where matter is highly clustered. Therefore, all studies 
(except for Ref. [13) use data on the largest scales. 

Using small scale data in galaxy surveys requires know- 
ing the nonlinear clustering and bias of galaxies relative 
to the dark matter. As we shall show, the precision of the 
small scale galaxy clustering data such as that from the 
SDSS is not high enough to warrant the inclusion of the 
effect of neutrino clustering in dark matter halos, though 
other systematic effects may be important (see Ref. JJjJ ) . 
However, when large weak lensing surveys - which mea- 
sure the mass distribution directly ~ become available, it 
will be essential to make direct use of this information, 
even on the smallest scales due to the expected precision 
of their results [l^ . 

The nonlinear power spectrum for the dominant clus- 
tering component, cold dark matter (CDM) itself has 
been best estimated by high resolution simulations by the 
Virgo Collaboration |16j . but their work quanitifies the 
uncertainty in their functional fits of the nonlinear power 
to approximately 7%. The effects of early free-streaming 
of neutrinos in suppressing the linear power spectrum 
can be incorporated into predicting the nonlinear power 
spectrum such as that from the fits of Refs. How- 
ever, one cannot naively expect this to characterize the 
full effects of massive neutrinos in the nonlinear regime. 

Vale & White 18] showed that the uncertainties intro- 
duced by approximations in ray-tracing techniques and 
numerical convergence of pure dark matter simulations 
may be sufficiently reduced with expected computing 
resources. Similar to the effects of neutrino infall into 
dark matter halos probed here. White ^^'^ Zhan & 
Knox showed that the effects arising due to baryonic 
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cooling and heating in CDM halos can alter the nonhnear 
matter power spectrum to significantly alter the observed 
weak lensing signal. Therefore, in order to eS^ectively use 
the information gained in upcoming weak surveys, one 
has to accurately determine the nonlinear matter power 
spectrum for a given cosmological model, and how it is 
influenced by the presence of baryons as well as massive 
neutrinos. 

In this paper, we first describe an analytic Boltzmann 
solution of neutrino infall into cold dark matter (CDM) 
halos to calculate the modification of captured neutrinos 
on the halos in We then employ the halo model to 
calculate matter clustering statistics including the effects 
of neutrino clustering, as well as the modification of the 
weak lensing power spectrum while including or ignor- 
ing this effect in ijllll and then sum up our conclusions 
in i jlVI In a companion paper |2H . we use multiparticle 
numerical simulations to quantify the effects of massive 
neutrino collapse into CDM halos on the nonlinear mat- 
ter power spectrum. 

II. CLUSTERING OF MASSIVE NEUTRINOS 
IN CDM HALOS 

To begin, we solve an isolated problem: how neutri- 
nos cluster in the presence of a dark matter halo. As 
we will see in the next section, the neutrino clustering 



around halos at late times leads to changes in the non- 
linear power spectrum. These changes are unique signa- 
tures of massive neutrinos and may eventually be de- 
tected. A pioneering work treated the clustering 
of massive neutrinos around a point-like seed with the 
Boltzmann equation. While this seed was taken to be a 
cosmic string, their technique was extended to accretion 
of neutrinos onto a CDM halo 23] . The density profile of 
a CDM halo has been found on average to follow a uni- 
versal profile over a wide range of mass scales, which we 
take to be a Navarro- Frenk- White (NFW) form [2ll25l| . 
The structure of the inner portion of the profile is not 
crucial to the neutrino clustering studied here. One as- 
sumption used in Ref. that simplifies the Boltzmann 
solution to this problem is that the NFW CDM profile 
is not influenced by the accretion of neutrinos; neutrinos 
do not act back on the CDM. (Formally, the source term 
in the Boltzmann equation is a pure NFW-type CDM 
distribution plus the neutrinos). This approximation is 
reasonable because the neutrino mass associated with a 
cluster is a small fraction of the mass in the cluster (of 
order 1%, see below) and more diffusely distributed, so 
that changes to the CDM NFW halo are essentially neg- 
ligible. 

In the notation of Ref. [22|, the Boltzmann solution 
for the Fourier transform of the neutrino profile in an 
evolving CDM halo is: 



Mk,d,,df) = "/3°;. , / dd'a\i)'){pcBM{k,^\MciyM{^')) + Mk,^^,^'})I{ ^^^^^^^^) , (2.1) 
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where a is the scale factor, but we parameterize integrals 
by a lookback super-conformal time, M = drj/a = dt/a?'] 
T^fi is the current neutrino temperature, 1.95K; CDM 
evolution is included in the term Mcdm ('?'); ^-ii*! 

//3 = / dx 2.2 
Jo -f 1 

The lower and upper limits of the integration set the 
time when the accretion starts, i^i, and when it stops, 
'df. The final time is chosen to be at the redshift that 
the power spectrum is needed. The initial time is chosen 
to be early enough such that an extremely high fraction 
of the neutrinos have velocities that cannot be captured 
by CDM halos, and deep enough to be beyond the galaxy 
distrubution of upcoming weak lensing surveys, at z = 
5. Our results arc essentially unchanged by choosing an 
initial z — 3. 

Before writing down the solution for the neutrino 
profile, we need to add one ingredient to the earlier 
work [23, 123 ■ We will be especially interested in the fea- 
tures of the power spectrum and its suppression around 



k — 0.5/1"^ Mpc, near the transition from the 1-halo to 
2-halo terms in the halo model (see below). Around 
k = 0.5h~^ Mpc, halos at and above 10^'^h~^ Mq dom- 
inate the power spectrum . These massive halos have 
collapsed only recently and were significantly growing 
during the accretion evolution history of the neutrinos. 
To include this, we need to insert a growth factor in the 
source term of the Boltzmann equation. The merging 
and grow th of CDM clusters was studied extensively in 
Ref. [23, who found that they evolve as: 

Mcdm{z) = AfcDM,todaye-^°'=^ (2.3) 

where the free parameter that defines the growth (ac) 
is given by the phenomenological relation M^,{ac) — 
0.018A/cDM, today HI] > whcrc is the characteris- 

tic nonlinear mass, defined where the fiuctuation scale 
^{■m,z) — 1 (see below). That is, Oc is the scale factor 
at which a halo of mass O.OISMcdm, today has collapsed. 
For very large k, on the other hand, very light clusters 
dominate the power spectrum and these clusters have 
collapsed in the distant past, so change very little during 
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the period from 2 = 1 until today, when most neutrino 
accretion takes place |2^. For increasingly large scales, 
the evolution of the CDM is significant: at the level of 
2%, 6%, 15%, 26% for lO^^, IQi^, lO", 10^^ solar mass ha- 
los, respectively. 

As it stands, the equation for the neutrino density 
transform H2.1|l cannot be simply integrated because the 
integral contains p^. This term in the integral equation 
represents the clustered neutrinos acting back on them- 
selves: as more neutrinos accrete onto a cluster, the clus- 
ter mass increases, and with it the ability to accrete more 
neutrinos. Since neutrinos make up a small fraction of 
the mass in a cluster, the effect of neutrinos pulling in 
more neutrinos can also be ignored as a second-order ef- 
fect. 

For cosmologies with more massive neutrinos that can 
make up a significant part of the cluster mass, the non- 
linear effects of the neutrinos acting back on themselves 
and on the NFW-type CDM cluster may become impor- 
tant. Ringwald & Wong [g^ found that the effect of 
this neutrino gravitational feedback is significant and the 
linearized Boltzmann equation approach underestimates 
neutrino infall by a small amount for less massive halos, 
but upto a factor of several for ~ lO^^M© halos, which are 
however rare. The magnitude of the effect as estimated 
here is therefore a minimum estimate of neutrino infall. 
The magnitude of the effect of infall will be quantified 
in detail beyond such 1-halo approaches in the multiple 
particle simulations in Ref. ^l"!. 

The inner /(/?) is sometimes approximated by letting 
the Fermi-Dirac denominator become a Boltzmann expo- 
nential 22]; then the integral can be performed analyti- 
cally. But Ref. 23] points out that this approximation is 
off by « 20% for large k, so we use the full Fermi-Dirac 
type distribution. 

The neutrino profile around a lO^'^h^^ Mq halo is 
shown in Fig. ^ We see that neutrinos do not cluster 
in the halos as efficiently as cold dark matter, so the 
profile inwards of 100 kpc is mostly flat and saturated 
(recall that the dark matter NFW profile increases as 
towards the center). The CDM profile is truncated 
at one virial radius. The choice to cut off the CDM at its 
virial radius suppresses the resultant neutrino population 
around 1 Mpc from the halo's center. 
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FIG. 1; A comparison of the spatial distribution of accreted 
neutrinos for an NFW CDM profile for a 10"/i"^Mq halo 
that has evolved in the fit of Ref. (2^ (solid line) and with a 
static profile (dashed line). The virial radius for this halo is 
4.8 Mpc for a an ^ = 0.7, Sic = 0.26 cosmology. 



neutrinos. The linear power spectrum determines how 
halos are correlated with each other. This first effect 
then is felt in the power spectrum in the so-called two- 
halo term, the contribution of halo-halo correlations to 
the total power spectrum. This first effect is due to the 
free-streaming of neutrinos on large scales, in the linear 
regime. After CDM halos form at redshift z < 5, massive 
neutrinos do not cluster as efficiently as cold dark mat- 
ter (again because of their velocities), but do fall into the 
potential wells of CDM halos. This second effect alters 
the profile of the matter within a given halo. 

A. Overview 

Let us first review the halo model [2^ |2^ . (See 
Ref. [33 for a recent review.) The nonlinear power spec- 
trum gets contributions from one- and two-halo terms: 
PNhik) = Pih{k) -t- P2h{k).', we first write them down 
and then explain the functions needed to compute them: 



III. THE HALO MODEL 

Here we show how to include the effects of massive neu- 
trinos within the halo model of nonlinear matter cluster- 
ing and then estimate their magnitude. The halo model 
posits that all matter exists in halos and the correlations 
in the matter can be explained by considering the cor- 
relations of the halos and the density profile of matter 
within the halos. In the context of this model, there are 
two places where neutrino masses affect the total mass 
distribution. First, the linear power spectrum in a model 
with massive neutrinos differs from one with massless 



and 

Pi.(fc) = J (3-2) 

Here p = Pcrit^c is the average matter density in col- 
lapsed halos the universe, i.e. flc is the fraction of the 
critical density that is in halos, including neutrinos that 
have fallen into halos. Note this is different from the 
usual halo model definition of p as the total average mass 
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density Pciit^m, where neutrinos are either not massive 
or ignored, and fl^ = ficDM + ^^6- Instead, we use 
flc = f^cDM + + ^^i^haio- ^uhaio is Calculated by in- 
tegrating the mass of neutrinos in individual halos over 
the halo mass function. 

The halo model integrates over regions with overden- 
sities parameterized by^ fj, = S'^/a{M)-^ , with 6c = 1-68 
the linear overdensity at the epoch of a halo's collapse. 
Here we neglect the effects of neutrino clustering on the 
definition of 6c- The rms of the linear fluctuations, <t{M), 
filtered on a scale which on average contains mass M is 



(M) 



(2^ 



PiUk)\WR{k)\' 



(3.3) 



where WR{k) is the fourier transform of a top-hat win- 
dow function of size R, the radius enclosing mass M. 
Each mass M then is associated with a particular value 
of II. Rare overdensities with large /x (high-cr peaks) cor- 
respond to large mass halos. The correspondence can be 
inverted to obtain M(/i), needed to compute the one and 
two- halo terms in Eqs. (|3.1|l and (|3.2|) . 

The number density of halos with mass M is deter- 
mined solely by the dimensionless ratio /i, which quanti- 
fies how rare the overdensity is; we use the Sheth-Tormen 
distribution Isal 



(3.4) 



where 



fit,) = A-{1 + {aii)-P)^exp{-afi/2), (3.5) 
M 

with A such that /(/i) integrated over all /i must be 
equal to 1 by mass conservation. We adopt a = 0.707 
and p = 0.3. The correlation between two halos de- 
pends on their mass. For cluster masses greater that 
M{fi = 1) = M*, halos will be much more strongly 
clumped than surrounding matter. This is the nonlin- 
ear clustering regime, whose threshold is indirectly set by 
the linear power spectrum. The halo model bias, 6(/x), 
folds in this mass-dependent correlation into the overall 
correlation of two halos, the "two halo" term. The bias 
associated with the Sheth-Thormen distribution is 



bifi) = 1 



2p 



^-1 + . 

6c 5c{l + {an)P) 



(3.6) 



B. Modification for Neutrino Clustering 

The components of the halo model described above 
can be readily updated if neutrinos have mass. They 



^ This ratio is usually denoted u, but we use fi here to save u for 
neutrinos. 
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FIG. 2: A comparison of the Fourier transform of the accreted 
neutrino profile into an NFW CDM halo of lO^/i'^M© that 
has evolved with redshift (lower line) and with a static profile 
(upper line). The Fourier transform for the CDM NFW profile 
is shown for reference of scale. 



depend solely on the linear power spectrum via cr, and 
one can easily insert the effect of massive neutrinos on 
the linear power spectrum. The final ingredient of the 
halo model is the mass profile around halos, represented 
by its Fourier transform p(fc,/i) in Eqs. (|3.1() and (|3.2(l . 
Usually, a natural choice is the NFW form mentioned 
above, 



r{r + r^)^ 



(3.7) 



Here r is the distance from the center of the halo (the 
conjugate variable to fc), M is the halo mass, and the 
two parameters ps and are functions of M and the 
concentration c. In terms of these, 



A„c3(M)p 



3[ln(l-fc(Af))-^^] 

3M 
47rc(M)3A„p' 



(3.E 



(3.9) 



and Ay is the virial overdensity with respect to the mean 
matter density, 



IStt^ + 82x - 39x^ 
1 + x ' 



(3.10) 



and X = flc{z) — 1 33j. We allow the concentration to 
vary with redshift z and halo mass 



c(M, z) 



9 



1 + z 



M 
If, 



-0.13 



(3.11) 



where M, is the cosmology-dependent characteristic non- 
linear mass scale. Cluster evolution becomes nonlinear 
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for virial masses greater than M*. The halo model is 
simplest if we work in Fourier space, where a given spher- 
ically symmetric density p{r) becomes 



P{k) = 



{inr'^ dr) p{r) 



sin(fcr) 
kr 



(3.12) 



We choose this cutoff to be at the virial radius of a given 
cluster Tvir = erg. 

While the NFW profile accurately describes the matter 
distribution if all matter is cold, it does not account for 
infall clustering of neutrinos. Since neutrinos have non- 
zero thermal velocities even at latest times, they will not 
strictly follow the cold dark matter profile in the halos. 
To account for this aspect of neutrino infall clustering, 
we need to generalize Eos. 13. II and 13. 21 by letting 



p{k,p) 
M{p) 



PCDuik, p) 



M,{p). 



(3.13) 
(3.14) 



the factor of 1 -I- z entering because r in Eq. (|3.12() is 
comoving distance. 

To include the effect of neutrino clustering, we add 
the clustered neutrino mass profile p^{k,p) to the halo 
model's fc-space density. In the spherical halo model, 
Pu{k,p) can be approximated by solutions to the Boltz- 
mann equation for neutrinos around an NFW CDM halo 
of mass M{p). Fig. |21 shows a comparison of the Fourier 
transforms of the NFW CDM halo profile and Boltzmann 
derived neutrino profiles. The fc — > limit of the Boltz- 
mann Eq. (|2.1(l gives the total neutrino mass accreted 
onto a cluster, M^{p), which appears in the halo model 
denominators, Eqs. (|3.17|I - H3.18|I . We do not include any 
potential changes the the form of the mass function f{p), 
since we assume the evolution of the CDM halos are un- 
perturbed by the presence of massive neutrinos. 

With these modifications, the nonlinear power spec- 
trum becomes 



Note that, for both cold matter and for neutrinos, 
limfe^o Pi{k) 



M,. 



(1 + 2)3 



(3.15) 



p^T^+^(fc) = p,^ 



CDIVH 
h 



'{k) + p: 



,CDM+iy 
2h 



(k), (3J6) 



and 

- P£r-W (/'"^ MJi:iTM,M fco„(t,rt^-P.(^.M)l)^ (3.18) 



We calculate the change in the nonlinear power spec- 
trum due to these modifications from neutrino clustering 
for a ACDM universe with parameters r^cdm — 0.26 — rijy, 
h = 0.7, fife — 0.04, n = 1 and erg = 0.9. Massive neutrino 
models are chosen with three degenerate mass neutrinos 
with nil/ = 0.1,0.3,0.6,0.9 eV (i.e., sum of all neutrino 
masses of 0.3, 0.9, 1.8, 2.7 eV), with chosen appropri- 
ately for these masses, while cts is fixed. The modification 
5PMk) = {P^i^{k)-P^^{k))/P^^{k) is shown in Fig.El 
Pj^L excludes the effects of late neutrino clustering and 
P]^L includes them. Both P^i^{k) and P^^{k) include 
the linear effects of early neutrino free streaming, since 
we are interested in the bias imposed by ignoring the 
effects of neutrino infall into CDM halos. The drop in 
power at fc ~ 0.5/i Mpc~^ occurs as expected at the scale 
where the most massive clusters are contributing to the 
nonlinear power spectrum, and increases with increasing 
neutrino mass. The reduction in power is due to the 
smooth structure of the accreted neutrino halo relative 
to the CDM halo. 



C. Weak Lensing Convergence Power Spectrum 

Planned weak lensing surveys have the potential to 
measure the power spectrum very accurately. There- 
fore, here we calculate the deviations due to massive neu- 
trino clustering on a weak lensing observable, namely the 
convergence power spectrum, Ci. This quantity is effec- 
tively the projected angular matter-matter power spec- 
trum weighted by the distribution of lensed galaxies. The 
signal is estimated by [sli |3^ [s^ 




(3.19) 

for a universe with flat geometry, where x is the comoving 
radial distance, Xh is the distance to the horizon, a = 
1/(1 -I- z), and P{k,z) is the nonlinear power spectrum 
at the appropriate redshift. The weak lensing weighting 
function is 

roc I 

9ix)^X dx'nix')^^, (3.20) 
Jx X 
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FIG. 3: The estimated change in the nonlinear matter power 
spectrum at z=0 calculated with the halo model, 5Pt^i^{k) — 
{PNh{k) — PSthik)) / P^i^{k) due to late time neutrino cluster- 
ing only. (Both P^i^{k) and Pt^i^{k), are calculated with a 
transfer function that includes the same early neutrino free 
streaming.) The lines of increasing magnitudes denote mod- 
els with three neutrinos with masses of 0.1 eV, 0.3 eV, 0.6 eV, 
0.9 eV, respectively. 



where n{x) is the redshift distribution of the lensed galax- 
ies normalized such that J dz n{z) = 1. 

The expected error in the observed weak lensing con- 
vergence power spectrum comes from two sources: sam- 
ple variance on large scales due to finite sky coverage, and 
on small scales by the finite number density of galaxies 
on the sky, 

= \/ r9f f^^ + ^) ' (3-21) 

V (2^+ l)/sky V "gal / 

where the fraction of sky covered /sky, intrinsic inferred 
galaxy ellipticity 7rms, and average number density of 
galaxies figai are survey dependent. 

As an example, we use survey parameters similar to 
those possible with the Large Synoptic Survey Telescope 
(LSST).^ Specifically, the galaxy redshift distribution 

is of the form n{z) cx z^e"^^/^"-* , with zo = 1, av- 
erage galaxy density ngai — 50 arcmin^^, a sky cov- 
erage of /sky = 0.5, and 7rms = 0.15. In Fig. ^ we 
show the effect on the weak convergence power spectrum 
SCe = {C^ - C^)/C^, when we include {C^) and exclude 
(C") the effect of neutrino clustering in CDM halos to 
see the dependence of the weak lensing signal. Note that 
the linear effect of early neutrino free streaming is in- 
cluded in both and C°, since we are interested in 
the bias imposed by ignoring the effects of neutrino infall 
into CDM halos. The cosmological parameters are cho- 
sen as in illll Bl Here, we take the neutrino infall effect 



^ http://www.lsst.org/ 
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2 

FIG. 4: The weak lensing convergence power spectrum (up- 
per panel) for nonzero neutrino mass models of 0.1 eV, 0.3 eV, 
0.6 eV, and 0.9 eV, with decreasing peak convergence, respec- 
tively. The power spectra are normalized at us — 0.9, there- 
fore showing a pivot at ^ ~ 200. The deviations including 
and excluding this effect are plotted in the lower panel, with 
increasing mass neutrinos corresponding to an increased am- 
plitude of the effect. Gray (cyan) boxes are expected errors 
for an LSST-like survey as described in the text. 



evaluated near the peak of the weak lensing weighting 
function g{x) at z = 0.4, causing a change in the shape 
of the weak convergence power spectrum deviation rela- 
tive to the z = nonlinear power spectrum above. The 
change in the nonlinear power spectrum affects the weak 
convergence power spectrum from ~0.1% for three 0.1 eV 
neutrinos to for three 0.9 eV neutrinos. 



IV. CONCLUSIONS 

We have estimated the effect of massive neutrinos in 
a concordance ACDM cosmology on the nonlinear power 
spectrum in the halo model. The shape of the nonlinear 
power spectrum is changed due to neutrino infall in CDM 
halos at a level of ~ 0.5% for three 0.1 eV neutrinos to ~ 
3% for three 0.9 eV neutrinos, corresponding to a change 
in the expected shape of the weak convergence power 
spectrum at a level of ~ 0.1% for three 0.1 eV neutrinos 
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to ~ 1% for three 0.9 eV neutrinos, being reduced due to 
the weight of higher rcdshift structures in cosmic shear 
measurements. 

Using the hnear information from weak lensing surveys 
we may be able to constrain the neutrino mass so that 
the effect of neutrino clustering may be neglected [38l |. 
Even given the inferred current upper limits on neutrino 
masses from the linear power spectrum (Stoj/. < 0.42 eV, 
95% C.L.) [13, the effect of neutrino clustering may be 
negligible unless current neutrino mass upper limits from 
the linear power spectrum are too stringent, or new de- 
generacies may emerge among features in the primordial 
power spectrum (e.g., running of the primordial spectral 
index or a break in the spectrum which allow 
for larger neutrino masses. 

For more general cases where massive neutrinos are 
present in fits to observed weak lensing convergence 
power spectra, massive neutrinos will need to be included 
in numerical simulation predictions of the weak lensing 
signal in addition to phenomena arising from baryonic 



condensation and heating 0, , leading to the poten- 
tial necessity of high-resolution multiparticle (neutrino, 
CDM and baryon) hydrodynamic numerical simulations. 
Coupled with upcoming weak lensing surveys, these pre- 
dictions will be a power probes of the contents of the 
cosmological soup as well as the process of cosmological 
structure formation. 
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